Canonical and non-canonical equilibrium distribution 
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We address the problem of the dynamical foundation of 
non-canonical equilibrium. We consider, as a source of di- 
vergence from ordinary statistical mechanics, the breakdown 
of the condition of time scale separation between microscopic 
and macroscopic dynamics. We show that this breakdown 
has the effect of producing a significant deviation from the 
canonical prescription. We also show that, while the canoni- 
cal equilibrium can be reached with no apparent dependence 
on dynamics, the specific form of non-canonical equilibrium 
is, in fact, determined by dynamics. We consider the special 
case where the thermal reservoir driving the system of interest 
to equilibrium is a generator of intermittent fluctuations. We 
assess the form of the non-canonical equilibrium reached by 
the system in this case. Using both theoretical and numerical 
arguments we demonstrate that Levy statistics are the best 
description of the dynamics and that the Levy distribution is 
the correct basin of attraction. We also show that the cor- 
rect path to non-canonical equilibrium by means of strictly 
thermodynamic arguments has not yet been found, and that 
further research has to be done to establish a connection be- 
tween dynamics and thermodynamics. 



I. INTRODUCTION 

A growing interest in non-extensive thermodynamics 
can be traced back to the 1988 pioneering work of Tsal- 
lis [Q. In particular, non-extensive thermodynamics is 
the proper theoretical context for developing the foun- 
dations of non-canonical distribution functions and non- 
canonical equilibria in physical systems. On the basis 
of the subsequent development of the original work by 
Tsallis, it has been argued that canonical equilibria, the 
distributions that form the basis of equilibrium statistical 
mechanics, are not generic. Rather, canonical equilibria 
are singular in a more general form of equilibrium, called 
a generalized canonical equilibrium |2|.[i| . A number of 
recent papers elaborate on this idea |4|-|6|], but we find ||] 
to be of special interest addressing as it does the founda- 
tions of both canonical and non-canonical equilibria. Ra- 
jagopal and Abe Q prove that the equilibrium described 
by the canonical distribution is not uniquely determined 
by the micro-canonical distribution, as one finds in text 
books. In fact if the phase space has a fractal, rather 



than smooth, structure, a non-canonical distribution will 
result. However, although canonical equilibria are not 
unique, and non-canonical equilibria are possible, one 
might conclude on the basis of the arguments in || that 
the form of the non-canonical distribution is uniquely of 
the form established by Tsallis and his co-workers jiJD. 
As attractive as the probabilistic and entropic arguments 
of Ref. p| are, they are not indisputable, and in fact we 
find herein, using dynamical arguments, that this is not 
the case. 

We show that the adoption of a dynamical approach 
to thermodynamical equilibrium yields a different con- 
clusion. First of all, we argue that the non-extensive 
condition based on memory, and probably that resting 
on long-range correlations as well, has the striking ef- 
fect of making the role of dynamics much more impor- 
tant than in the case of ordinary statistics. The dynamic 
approach, in these non-extensive conditions, generates 
a form on non-canonical equilibrium that, however, de- 
parts from the generalized canonical form prescribed by 
non-extensive thermodynamics. To make it easier for 
the reader to follow our arguments and to und erst and 
the purpose of the paper, illustrated in Section ID, we 



shall first discuss the different natures of the entropic, 
dynamic and stochastic approaches to equilibrium. 



A. Non-extensive entropic indicator 

The argument for the non-extensive thermodynamics 
of Tsallis goes as follows. First of all, the conventional 
entropy of Gibbs 



S(U) = - J dx II(a:) lnU(x), 
is replaced by the non-extensive entropic indicator 

i-n(x) 9 



S q (U) = J dr 



1 



(1) 



(2) 



Secondly, we have to apply a method of entropy max- 
imization under given physical constraints to determine 
the most plausible shape of the unknown probability den- 
sity function II(a;) Q|. The first constraint is on the nor- 
malization of the distribution II(x): 



dx U(x) = 1. 



(3) 
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The second constraint is on the first moment of the vari- 
able x itself. According to the most recent prescription 
of Ref. H the constraints on x must be applied ||] on the 
mean value U q defined by 



fdxx U(x) q 
U q = ± 



/ dx Il(x)i 



(4) 



It has to be pointed out that in Ref. |g] the physical mean- 
ing of x is that of energy. Here, we shall interpret x as 
the "coordinate" of an overdamped particle driven by a 
fluctuation-dissipation process resulting from the inter- 
action with a non-conventional "thermal bath" . When 
the friction term can be neglected we are expected to re- 
cover the results of the earlier work of Ref. |1(J . Thus, for 
the same reasons as those illustrated in Ref. jlO| we 
set a constraint on the first moment of \x\. The result of 
entropy maximization subject to the imposed constraints 
of Eq. (|) yields 



n(x) = 



1 - 



(1 - q)b(x - U q ) 



Jdx n(x)« 



1/(1-3) 



/z q 



(•5) 



where Z q is a normalization factor and (3 is a constant 
value stemming from the Lagrange multiplier associated 
with the constraint on the variable x. Note that the ex- 
pression of n(x) provided by the authors of Ref. || is not 
explicit. In fact, as shown by Eq.(||), it is a functional of 
II(a;). Thus, in practice the explicit form of H(x) has to 
be established by means of an iteration procedure. It is 
remarkable, however, that the micro-canonical derivation 
from Ref. Q results in the same prescription as that of 
Eq. ([5]) . From the point of view of the issues under dis- 
cussion here what matters is the fact that at the end of 
the iterative procedure U q becomes a well-defined num- 
ber. Consequently, the resulting expression for H(x) is 
a simple analytical formula that in the asymptotic limit, 
|x| — > 00, with the constraint on the first moment of \x\, 
has the same structure as that derived in Ref. (lOl, 



11(1) 



6(2 - q) 



[l + 6(g-l)|ar|]V(9-i) 



(6) 



Note that the adoption of a constraint on the second 
moment would lead to 



U(x) 



[5(9-1)" 


r(l/(g-l)) 


7T 


r((3- 9 )/[2(g-l)j) 



(7) 



Note that both Eq.(^) and Eq.(0) are the result of 
an earlier prescription, missing the normalization factor 
present in Eq.(4). Nevertheless, it has to be pointed out 
that both Eq.(6) and Eq.(||) share the characteristic of 
having long tails with Levy statistics. In fact, it is well 
known that the anomalous diffusion processes of the 
Levy kind, in the one-dimensional case, are characterized 
by probability distributions p(x, t) whose Fourier trans- 
form in the symmetric case reads 



p(k,t) =exp(-b\k\ a t), 



(8) 



where a is the Levy index ranging, in principle, in the 
interval < a < 2, and b denotes the diffusion inten- 
sity. The inverse Fourier transform of p(k, t) of Eq.(||) is 
characterized by the tail |l| 



lim p(x, t) cx 



\l+a 



(9) 



which would lead immediately to the anomalous entropy 
index 



l + l/(l + o). 



(10) 




FIG. 1. Comparison of a Levy distribution (solid line) 
with distributions obtained by maximizing Tsallis entropy 
(dashed lines). The Levy curve is obtained by the inverse 
Fourier transform Levy characteristic function with a = 1.5, 
b — 4.4 • 10 4 . The asymptote is proportional to 1/x 2 ' 5 . The 
long-dashed line comes from Eq. (pi) with the constraint on 
the first moment (b = 0.004, q = 1.4). The short-dashed 
line comes from Eq. (pj) with the constraint on the second 
moment (b = 1.4 • 10~ , q = 1.8). The values of the param- 
eter b are selected, in both cases, so as to fit the asymptotic 
behavior of the Levy distribution. This constraint results in 
significant departures from the Levy distribution at small and 
intermediate distances. 

We have to stress, however, that except for the case 
a = 1, corresponding to the ballistic limit, the central 
part of the distribution produced by the Levy process 
of Eq. (|^) cannot be expressed in an analytical form and 
can smnificantly depart from the analytical form of both 
Eqs.jH) and (J5[) . The adoption of the prescription of 
Eq.([7f) is expected to yield better agreement, but, as we 
shall see subsequently, it does not seem to be compatible 
with the nature of the dynamical approach to equilib- 
rium. In Fig.|l] we compare the Levy distribution to both 
the prediction of Eq. (^) and Eq. (Q) and find that there ex- 
ists a significant disagreement between Tsallis and Levy 
statistics, even though some authors || refer to them as 
equivalent. In this paper we focus our attention on the 
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origin of the difference between the two kinds of statis- 
tics. Our dynamic approach leads to a form of equilib- 
rium that is stable, in the sense of the Levy-Gnedenko 
theorem |3 , while the generalized canonical equilibrium 
of Eq.(H) is not. This is evident in the free case, due to 
the difference between the Tsallis and the Levy structure. 
In Section |l| we prove that this is so also in the presence 
of a feedback with the system of interest of the generator 
of fluctuations. 



B. From dynamics to thermodynamics 

We are convinced that there are no incontrovertible 
reasons why the canonical distribution should be the 
unique form of thermodynamic equilibrium, and on this 
issue we essentially agree with the point of view of Ra- 
jagopal and Abe However, we are equally convinced 
that the generalized canonical form of Eq.fg) does not 
satisfy the stable conditions necessary for Levy statis- 
tics, and the results of the present paper can be thought 
of as providing plausible evidence of that. To substanti- 
ate this view with dynamical arguments, it is convenient 
to concisely review the results of an earlier work 

The ambitious purpose of this earlier work was that of 
reversing the path from thermodynamics to mechanics 
established by Boltzmann. The main idea behind Ref. 
|l5| is as follows. The Fokker-Planck equation is a well- 
known description of the evolution of the probability den- 
sity in phase space. Many attempts have been made in 
the literature to derive this equation Jl6| pc|. However, 
all these attempts rest on the assumption that the bath, 
responsible for the Brownian behavior of the particle of 
interest, is given by a set of harmonic oscillators. This 
means that a statistical assumption must be made on the 
initial condition of the bath, which is arbitrarily given 
a canonical equilibrium distribution, corresponding to a 
given temperature T of the thermal bath. 

Consequently, this kind of approach to Brownian mo- 
tion is only partially dynamical, since significant use is 
already made of statistical mechanics, and thus of ther- 
modynamics The authors of Ref. |l5| adopted a 
totally different approach. They assumed that a given 
oscillator, playing the role of a stochastic Brownian par- 
ticle, interacts with another Hamiltonian system, which 
should play the role of a bath. They called this second 
system a booster, to stress that the ensuing approach has 
to rest only on the dynamical properties of this kind of 
bath, with no use whatsoever of thermodynamical argu- 
ments. After establishing the dynamical conditions en- 
suring the validity of the Fokker-Planck equation, it is 
found that the oscillator of interest reaches a canonical 
equilibrium distribution. Due to the nature of the pro- 
cedure adopted the width of this canonical distribution 
depends only on the parameters of the Hamiltonian sys- 
tem under study. Consequently, it is possible to derive a 
mechanical expression for temperature. This is the key 



result of Ref. |15|, which reads 
k B T 



^lnA(E) + JLln{(e)e q Re[^(w)}} 



(11) 



Note that £ denotes the doorway variable, namely the 
variable of the booster through which the interaction be- 
tween the booster and the oscillator of interest is estab- 
lished. The symbol denotes the Laplace transform 
of the correlation function of £ evaluated at the oscilla- 
tion frequency of the oscillator. The structure of this ex- 
pression reflects the application of linear response theory 
pq| . The correlation function whose Laplace transform 
is in Eq. ( |il"| ) is evaluated assuming the booster to be 
in a micro-canonical equilibrium with energy E and this 
condition is not affected by the interaction with the os- 
cillator. The symbol A(E) denotes the number of states 
of the booster in the same physical condition, and, con- 
sequently, obeys the ordinary prescription |22| 



A(E) oc E N ' 2 



(12) 



where N is the number of degrees of freedom in the 
booster. 

The authors of Ref. [|l5| note that for N — > oo, 

—ln{(e) &q Re[^)}} « gE lnA (E)- (13) 

This means that in the limiting case of infinitely many 
degrees of freedom the predictions of Boltzmann are re- 
covered: 1/(/cbT) = -^lnA(E). Note that the condition 
of Eq. ([l3]) holds true if the doorway variable £ depends 
on a number of particles which is kept fixed with N — > 00. 
In other words, the condition of Eq.([l3|), or equivalently, 
Boltzmann's principle, rests on an interaction condition, 
which is crucial for the extensive statistical mechanics 
perspective to apply. 

Let us see this aspect in detail. In Ref. Jlj| the booster 
is the the well known Fermi-Pasta-Ulam (FPU) system 
p3| . This is a chain of particles interacting with one an- 
other via nonlinear interactions. Let us consider two op- 
posite conditions. In the former, which is the one adopted 
in Ref. [p| , the oscillator of interest interacts only with 
the first particle of the FPU chain. In the latter the os- 
cillator of interest interacts with all the particles of the 
FPU chain, with interaction strength of comparable in- 
tensities. The former condition refers to a short-range 
interaction, confined to the positions of the oscillator of 
interest and to the first particle of the FPU chain. In the 
latter condition, the interaction extends over the whole 
FPU chain. It is evident that the former condition fits the 
inequality of Eq.(|T3|), whereas the latter does not. This 
means that the latter condition results in a non-extensive 
form of dynamics, with a consequent breakdown of the 
prescriptions of ordinary statistical mechanics. The for- 
mer condition, on the contrary, for N — > 00, recovers 
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ordinary statistical mechanics. This suggests that the dy- 
namical corrections to the Boltzmann principle, recorded 
in Ref. Gsl for relatively small values of N, are a mani- 
festation of incipient non-cxtensive statistical mechanics, 
and so are very close to the breakdown of the Fokker- 
Planck treatment on which the analysis Ref. rests. 
Note that when we assign to the oscillator of interest a 
very low frequency, compared to the booster frequencies, 
the quantity Re[^^(u)] turns out to virtually coincide 
with the time scale T m i cro of the variable £, defined by 



/>oo 

T m icro = / 3>f(£) dt < OO. 

Jo 

Herein we focus our attention on the case where 



with 



lim t —, 00 '&( i {t) = const/r, 



<0 < 1. 



(14) 



(15) 



(16) 



The inverse power law form of the correlation function 
means that we select time memory as the source of vi- 
olation of extensivity, rather than long-range spatial in- 
teractions. The parameter T m i cro denotes the correlation 
time of the fluctuating variable £. Here we consider the 
dichotomous case, where the variable £ has two distinct 
values, W, and —W, with fluctuating time durations. 
The correlation function $^(i) is proven [Q to be pro- 
portional to the second time derivative of the distribution 
of waiting times in the two states of the variable £. This 
function, as we shall see in Section j~C| depends on an- 
other time, T, as well as on j3. 

In conclusion, the work of Ref. Jl5| establishes the dy- 
namical conditions necessary to derive canonical equi- 
librium. The authors prove that canonical equilibrium 
implies an interaction with a booster with a finite time 
scale and a short-range interaction. We can also observe 
that in this case the use of dynamical arguments yields 
the same conclusions as the very simple argument based 
on the law of large numbers. The advocates of the law of 
l arg e numbers for the foundation of statistical mechanics 
|p5[ might judge the dynamical perspective to be of lim- 
ited use. We prefer to interpret the conclusion of Ref. 
Jl5| as proof that in the extensive case the canonical 
distribution description of thermodynamic equilibrium is 
unique. The numerical work of Ref. |ll| shows that it is 
very difficult to detect the dynamical corrections to the 
Boltzmann principle illustrated by Eq.(fil"|) from within 
an extensive perspective. This is so because the dynam- 
ical corrections become significant when the booster is 
small, and, as clearly pointed out by the recent work of 
Gross ^6|, a small system is essentially non-extensive. 
Thus, the numerical findings of Ref. [[15[ refer to a con- 
dition very close to the breakdown of the extensive con- 
dition and consequently of the canonical equilibrium on 
which Eq.(pl]) rests. This means that the dynamical ap- 
proach can be of great utility. In the case of non-extensive 



statistical mechanics, the dynamic approach is probably 
the only non-ambiguous way to address the subtle issues 
posed by the entropic and probabilistic methods. 



C. Stochastic dynamics 

In this section we explain the nature of our dynamical 
approach to the non-canonical equilibrium distribution. 
This approach rests on a stochastic method adapted to 
the need of realizing a dichotomous variable £ with the 
two possible values, W and — W , and with the waiting 
time distribution 



m = 



(/3T + i) 2+/5 ' 



(17) 



where ranges in the interval of Eq.(|T^). Note that 
the parameter T is of crucial importance to define the 
time scale of our process and corresponds to the mean 
residence time in either of the two states of the veloc- 
ity variable £. In accordance with the prescriptions of 
Ref. (24) this waiting-time distribution yields the kind of 
correlation function ^(i) that we plan to study herein 
(see Eq.(|l5|)). In fact, as shown in Ref. ]24j] , the form of 
Eq. ( |l7| ) yields the correlation function 



d> c (i) 



(pry 



{I3T + t)f> 



(18) 



fitting the asymptotic time limit of Eq.(|l5|). This con- 
dition can be considered as the natural one-dimensional 
counterpart of the two-dimensional billiards of Zaslavsky 
p7f . From this point of view, our commitment to the 
adoption of a merely dynamical approach is not broken, 
since the stochastic approach that will be adopted in Sec- 
tions [n] and [II is statistically equivalent to the adoption 
of the intermittent map of Ref. |2|] , on which the theo- 
retical work of Ref. |n| is based. 

We stress that the equivalence between a dynamical 
map and a stochastic process is the reason why contact 
can be established between the dynamical and the en- 
tropic approaches. In fact, as shown in Ref. (h]] the non- 
extensive Tsallis entropic indicator serves the purpose 
of guessing the most convenient form for the transition 
probability Tl(x), which is then related to the waiting- 
time distribution ip{t) through the basic property 



n{x) = ^ip(x/w)/w. 



(19) 



The factor of 1/2 takes into account that the probability 
of making the jump x is equal to that of making the jump 
—x, namely, the jump in the opposite direction. The au- 
thors of Ref. flTj|| show that the left-hand term of this 
equality, Eq.(|l9|), can be predicted using entropic argu- 
ments, while the right-hand term of the same equation is 
dictated by dynamical arguments based on the intermit- 
tent map of Ref. |2q]. These dynamical arguments are 
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supplemented by the assumption of random injection of 
the trajectory from the chaotic into the laminar region 
of the intermittent map, an argument leading to an an- 
alytical prediction for ip(t) in complete agreement with 
the numerical observation of dynamics p8j. 

The earlier work of Rcf. |ly] established that the adop- 
tion of the method of entropy maximization applied to 
the non-extensive entropy of Eq.(||) results in a form of 
H(x) which is compatible with the birth of Levy statis- 
tics. However, the H(x) thus derived is not the equi- 
librium distribution of the variable x. Rather it is the 
probability for the random walker to make a jump of 
length \x\. This is a stationary property determined by 
the special kind of booster here under study. The ensu- 
ing diffusion process yields a Levy form as a result of the 
Levy-Gnedenko theorem ]l4|] . As shown in Fig. [I], the 
shape of this distribution departs from the form of the 
generalized canonical distribution of Eq.(||), even if we 
adopt the constraint on the second moment, in spite of 
the fact that it does not fit the nature of the dynamical 
approach here illustrated. 

To account for this discrepancy we might make the 
conjecture that the comparison between the Levy statis- 
tics and the Tsallis generalized canonical distribution is 
not appropriate. The former refers to a diffusion process, 
and the latter to an allegedly equilibrium condition. Ac- 
tually, after exploring this possibility we shall conclude 
that the latter, at least in the case of the dynamic model 
of the present paper, cannot reflect an equilibrium prop- 
erty. However, at the present stage, we are forced to 
develop a picture comparable to that of the earlier work 
of Ref . [|l5| . We have to study a case where the variable x 
not only undergoes the influence of the diffusion produc- 
ing fluctuations, but it produces a feedback on its own 
"bath" , balancing the diffusion process, so as to create an 
equilibrium condition. This is the condition to compare 
to the generalized canonical distribution of Eq.(|^). In 
other words, to address, from a dynamical perspective, 
the issues recently dealt with by Rajagopal and Abe || 
we cannot disregard the feedback from the system to the 
"bath." 



D. Purpose of the present paper 

At this stage it is much easier for us to illustrate the 
main purpose of the present paper. We want to explore 
a condition where the booster does not fulfill the key 
condition of Eq.([l4|), that is, the microscopic relaxation 
time diverges, and we want to assess whether or not this 
divergence leads to the generalized canonical distribution 
of Eq.Q. We aim at answering the question: Which is 
the form of the equilibrium distribution reached by the 
system of interest when the time scale of the booster is 
infinite? 

To answer this question we do not adopt, as done in 
Rcf. a Hamiltonian approach. The latter approach 



is difficult for obvious reasons. A numerical simulation 
check of the theoretical prediction would imply technical 
difficulties caused by the slowness of the booster itself. 
This means that the anomalous booster is replaced, as 
done in earlier papers |29|]3C| ], by a generator of a di- 
chotomous fluctuation with an inverse power law distri- 
bution of waiting times. The essential ingredient of the 
approach of Ref. |fj"5|j , not yet present in the dynamical 
derivation of a free Levy diffusion, is a feedback of the 
diffusing variable on the generator of the fluctuations. 
This aspect has already been considered in the earlier 
work of Ref. p9|] . However, in that paper the dynamical 
approach to the feedback was replaced by a phenomcno- 
logical friction, which did not allow the authors of that 
paper to keep the promise of resting solely on dynami- 
cal arguments at any level. In conclusion, we adopt the 
program of Ref. JL5J, based on the observation of the 
fluctuation-dissipation process caused by the interaction 
of a particle with a booster having no finite time scale. 
The presence of feedback serves the purpose of balanc- 
ing the diffusion process with dissipation so as to result 
eventually in an equilibrium condition. 

With the program of Ref. Jl5| in mind, we have to re- 
fer ourselves to the correlation function, Eq. ([18]) , with the 
index f3 fulfilling the condition of Eq. ( plf ) , and so imply- 
ing the breakdown of the time scale separation between 
the macroscopic and the microscopic levels. The real- 
ization of the program of statistical mechanics requires 
an accurate definition of the process of memory erasure 
associated with the transition from one to the other ve- 
locity state This is more conveniently defined by 
the waiting-time distribution ip(t) of Eq. (|l7|) , than by the 
correlation function of Eq. (jl^) . We see that even if the 
condition of Eq.(|l6|) applies, the time f = t ?p{t)dt 
remains finite. Memory of microscopic dynamics is lost 
in times t >> T. 

At this stage it is convenient to support our claim 
concerning memory erasure with arguments borrowed 
from the earlier investigation of Gaspard and Wang [^2| . 
These authors studied the Kolmogorov complexity of the 
Manncvillc map and found that in the regime correspond- 
ing to the dynamical foundation of Levy processes the 
Kolmogorov complexity is a linear function of time. This 
means that at a given time t the number of transitions 
M from the one to the other laminar region is given by 
M (xt/T. This means that for t — > oo the conditions for 
the realization of the generalized version of the central 
limit theorem are fulfilled, since the function defined by 

p{x,M) =n(x) *n(a;) * . . . * U(x) *p in (x), (20) 

where the asterix denotes a convolution, and for M — > oo, 
Eq. ( pp[ ) tends to the Levy distribution. One might be 
tempted to make the conjecture that in the presence of 
feedback the number M cannot increase beyond some 
limit, and that equilibrium is reached with a relatively 
low value of M so as to allow the distribution to maintain 
the structure of a generalized canonical distribution d 
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la Tsallis. In this paper we limit our analysis to the 
case where the linear response theory of Section IE holds 
true. This forces us to adopt a feedback so weak that 
this possibility does not emerge from our simulations. 



II. DYNAMICAL MODEL 

The dynamical model studied here depends on a pro- 
cess of free diffusion with feedback, established through 
control of a dynamical parameter of the booster. This 
feedback process has the role of balancing diffusion so as 
to realize dynamically the equilibrium condition. We il- 
lustrate first the dynamical model used as the generator 
of free diffusion, then we show how the feedback is re- 
alized, and finally we illustrate the numerical technique 
adopted. 



A. Free diffusion 

Formally free diffusion is realized using the equation of 
motion 



dx(t) 
dt 



£(<)■ 



(21) 



Here the dynamical variable x(t) denotes either a spatial 
coordinate or a velocity. In the former case the variable 
£ (t) has to be considered a fluctuating velocity, while in 
the latter case it is regarded as a fluctuating acceleration. 
The results are equivalent, and the reader can adopt ei- 
ther of them, even if we consider x(t) to be a spatial co- 
ordinate so as to make the connection with earlier work 
JTo| more natural. 

We assume the variable £(t) to be dichotomous, 
namely, we shall assign to this variable only two distinct 
values, either W or — W. The motivation for this choice is 
not only simplicity, but also has to do with the main pur- 
pose of this paper, that of assessing the specific form the 
non-canonical equilibrium. Also, of course, and prior to 
this purpose, that of proving that deviations from canon- 
ical equilibrium can be generated as a consequence of the 
breakdown of the ordinary condition of time-scale sepa- 
ration between macroscopic and microscopic processes. 
In the limiting case where the variable £(t) is Gaussian, 
and consequently produces a Gaussian equilibrium dis- 
tribution p3| , the resulting diffusion process is Gaussian 
and the final equilibrium distribution can turn out to be 
Gaussian as well J2j| , as a result of having established a 
seed of ordinary statistics at the microscopic level. For 
this reason the choice of microscopic statistics, depart- 
ing from the ordinary statistical condition is crucial, and 
the adoption of the dichotomous assumption serves the 
purpose of eventually establishing a non-canonical equi- 
librium. 

The variable keeps one of the two possible values 
for times with a random duration. Thus, a statistical 



treatment must be adopted and a waiting-time distribu- 
tion ip(t) is used. The time intervals of sojourn in a given 
state are labelled by an integer index k running from 1 
to +00. At the end of any sojourn the variable can cither 
change its value, and thus make a transition from W to 
— W or from — W to W, or it can keep its original value. 
The probability of changing values and that of keeping 
the same value are equal, and thus equal to 1/2. We can 
interpret the variable £ (t) as a function of the continuous 
time t by adopting the following prescription: 



(22) 



where the times tk and tf-+i denote the beginning and the 
ending time of the k — th sojourn, respectively, and the 
values £fc are randomly assigned either the value W or 
the value — W with equal probability and Xt k ,t k+1 (t) = 1 
if tk < t < tk+i- The sequel of the times t' k is fixed by the 
waiting time density distribution ip, by selecting a given 
value Tk , of the interval [Tk , Tfe + e] , with e < < 1 with 
the probability ^(Tfe)e. The waiting time distribution ip 
is assigned an inverse power law form determined by the 
constraint of yielding the correlation function 



*«(*) 



(pry 



ipr + ty ' 



(23) 



where (3 is a positive number determining the integrabil- 
ity of the correlation function. If (3 > 1 the correlation 
function, Eq.([2^), is integrable and the microscopic time 
scale 



(3T 



10 P ~ 1 



(24) 



can be defined. According to Refs. |54|,[24j the form of 
ip(t) is determined by that of <&{(£) through the relation 



which yields 



(25) 



(26) 



We note that the parameter T appearing in both Eq. ( |25 
and Eq.(|26|) is the mean sojourn time. Thus, in some 
sense, when the microscopic time of Eq.(p4]) becomes in- 
finite, the role of the microscopic time scale is played by 
T. 

By integration of Eq. (Glh we get 



!(t) =&,(«- *n-l) + X) 6 



kTk, ■ 



(27) 



fe=0 



where t is a time located in the interval [i n _i, £ n -i + T n] 
with t„_i = yV,_ Q Tfe, in accordance with the earlier 



G 



prescrip tion s. Furthermore, the transition probability of 
Section IA is obviously related to the waiting time dis- 
tribution by Eg. (|I9|) . Note that, as discussed in detail in 
Refs. ^,^,^5|, there are two distinct basins of attrac- 
tion for the diffusion process resulting from the repeated 
occurrence of the transition process of Eq. , the Gauss 
basin for (3 > 1 and the Levy basin for (3 < 1. 



B. The booster linear response and the diffusion 
control 

To control diffusion in such a way so as to generate a 
known equilibrium distribution we adopt a model based 
on replacing Eq.(pll) with 



dx{t) 
dt 



(28) 



which means that the motion of the fluctuating variable 
S,x(t) also depends on the value of the dynamic variable x. 
We make the simplifying assumption that the relaxation 
time, or response time, of the booster, Tb, is shorter than 
that of the particle of interest, Tr, so as to ensure for 
the booster a time dependent equilibrium condition de- 
termined by the state of the system of interest. In other 
words, the booster is assumed to be in a condition of equi- 
librium determined by the variable x. As the variable x 
moves from an initial condition x(0) out of equilibrium 
state, i.e. from a position larger than that fluctuations of 
x itself, to the hnal equilibrium, the booster correspond- 
ingly moves from an equilibrium to another equilibrium 
condition. All this is illustrated by the following mathe- 
matical arguments. These arguments do not aim at pro- 
viding a rigorous mathematical treatment, but rather a 
heuristic treatment reflecting the physical condition that 
we are assigning to the booster. 

To make the analysis in terms of a strictly stochastic 
process, we build up the following regularized process 
XT p (t), defined by: 



dx Tp {t) _ x(t)-x(t-T p ) 



dt 



t, x {t')dt' 



P Jt-To 



^(t')dt'. 

1 P Jt~T„ 



(29) 



where T p is a regularization time and in the last equiva- 
lence it has been explicited the spatial dependence of the 
noise on x. This new process x has the characteristics of 
x(t) in the mean, but its fluctuations are smoothed by 
using the temporal mean. 

If the response time Tg of the bath is fast and if T p is 
long enough, but always smaller than relaxation time Tr 
of the variable x, one can suppose that the phase space 
mean is almost equivalent to the temporal mean: 



dx Tp (t) 
dt 



1 



Tn 



(30) 



Now, assuming linear response is valid, the mean value 
of the thermal bath variable depends on x as: 



(31) 



With the linear response hypothesis inserted into Eq. 
(p0|) the equation for the regularized process 



dx Tp jt) 
dt 



-lx Tp {i) 



(32) 



is obtained. The solution XT p {t) = x(0)e 7 * shows a 
relaxation toward the equilibrium x ~ 0, on a time scale 
Tfl~l/ 7 . 

The analogy with a stochastic equation, such as the 
Ornstein-Uhlenbeck or its extension to Levy noise , is 
completed by adding a fictitious stochastic noise d£,r p (i) 
in (|32|), that summarizes the fluctuations on the time 
scale T p and that compensates for the approximation in- 
serted into Eq. (|30"|): 



(33) 



dx Tp (t) = -iZ Tp (t)dt + d£ Tp (t) 



Hence, with the assumptions of linear response theory 
and relaxation times of the thermal bath faster than of 
that of the variable of interest (partial equilibrium), a 
diffusive process becomes an equilibrium process. We 
emphasize that the analogy (p3|) is true only if: 



T B «T p «T R = 1 - 



(34) 



By using this simplified analysis it is also possible to 
make qualitative predictions about the equilibrium dis- 
tribution. 

When x ~ 0, the feedback —jxt {t) can be neglected 
with respect to the fluctuations £r (t). This limit corre- 
sponds to free diffusion for which the central limit theo- 
rem (generalized or not) can be used. The distribution is 
Gaussian or Levy's according to whether the correlation 
function is integrable or not. 

When | a; | >> the feedback term ( |3l| ) prevails over 
fluctuations in Eq.(j33|). Then the dynamical variable is 
pushed back to values near the equilibrium x ~ as an 
attractive field and, so, quenching the diffusion. This 
description corresponds, at least for the region \x\ << 
W/j, to what one could obtain if in Eq. ( |33| ) the fictitious 
noise £t p (0 is simply stochastic p9[ . 

The condition ( pl| ) is satisfied when the time evolution 
of £ x is rep laced by the stochastic variable appearing 
in Eq. ( [22|) with the the stochastic variable £,k(xk), which 
has two distinct values, with the same probability: 



W+(x k ) = W- 7 x k 
W-(xk) = -W - 7x fe 



(35) 



In the time scale T p , we can replace the variable £ x (t) of 
Eq.(p8[) with its average over the faster bath fluctuations, 
and thus with 



t-T„ 
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jx — (W + jx)] = —jx. 



(36) 
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FIG. 2. Sample trajectory for the diffusion process with 
feedback defined in Eq. (|2£|). The solid and dashed lines 
denote the states W+ and W-, respectively. See Eq. (|3^) for 
the definition of these two states. The two horizontal lines 
defining the central stripe correspond to the levels W/'y and 
— W/f. The parameters used have the values: 7 = 5- 10~ 3 , 
W = 1,(3 = 0.5, T =50. 

In Fig. U we show a sample trajectory corresponding 
to the dynamic prescriptions of Eq. (^) . In the lateral 
region, after the end of the state that has pushed the tra- 
jectory outside the central region, the slope assigned to 
the trajectory by the ensuing states have always the same 
sign, either negative, for the top region, or positive, for 
the bottom region. This is a consequence of the definition 
itself of side regions and central stripe, the former being 
those characterized by |x| > W/j and the latter referring 
\x\ < W/7. The final effect is that the trajectories that 
have exited the central region, are steadily brought back 
to it. We see that, on the contrary, in the central region 
the trajectory fluctuates back and forth more or less as 
it would do in the case of free diffusion. This is per- 
haps the intuitive reason why the linear response theory 



of Section III, referring to a central region of extremely 
large size, makes the system's statistics fall in the Levy 
basin of attraction as it does in the free case |3Ql. 



C. Computer simulation 

The computer simulations are done generating the tra- 
jectories of Eq.(27) with the assumption that the values 
W + (x k ) = W — jxk and W-{xk) = —W — "fXk have the 
same probability. The distribution of sojourn times in 
one of these two states is given by the function ip(t) of 
Eq. (p6[) , which, in turn, is realized by a suitable nonlinear 
transformation of a variable with a uniform distribution 
in the interval [0, 1] obtained by a standard routine for 
the generation of random numbers. It has to be pointed 



out that the value of either W+(xk) or W-(xk) is deter- 
mined by the value that the variable x has at the moment 
when a new state is established and it is kept fixed until 
the end of the sojourn in this state, namely its value is 
fixed in the interval [tfc , ifc+i] ■ Each trajectory is gener- 
ated according to the integral of the equation (^) with 
the initial condition a;(0) = and it is observed until 
the time Tf = 2O/7 >> Tr. At this time the value is 
recorded in the form of a histogram. The distribution of 
the variable x at this final time, assumed to correspond to 
equilibrium, is obtained by generating a sufficiently large 
number of trajectories and constructing a histogram. 



III. THEORETICAL PREDICTION 

We are now in a position to demonstrate that in the 
case of extremely weak feedback the resulting equilibrium 
distribution is Levy. There are no compelling theoreti- 
cal predictions in the case of strong feedback, except the 
heuristic arguments used in section |. 

According to the strategy adopted in Ref. jl5| the the- 
oretical approach to friction implies the treatment of the 
response of the booster to an abruptly applied external 
perturbation. Using the notation adopted in Ref. jD| the 
solution to the equation of motion 



^,t) = T{K)p^,t), 



(37) 



where T(K) denotes the operator driving the motion of 
the probability distribution p(£, t) in the presence of a 
perturbation of intensity K. We assume that this oper- 
ator can be expressed as the sum of an unperturbed and 
a perturbed part, 



T(K)=T + KT 1 , 



(38) 



where Tq denotes the operator driving the unperturbed 
motion of the variable £ and KTi denotes the operator 
corresponding to the perturbation. The intensity of this 
perturbation can be freely changed by modifying the size 
of the parameter K. We apply the perturbation abruptly 
at time t = 0, to the booster, assumed to be in the equi- 
librium state, p eq . The perturbation strength is assumed 
to be so weak as to make it possible to solve Eq. (|37|) at 
the first order perturbation, 



^pi(U) =T p 1 (£,t)+KT lPo , 



(39) 



where po and p\ denote the zero-th and the first-order 
distribution density, respectively. We assume that the 
system before the abrupt perturbation is in thermody- 
namic equilibrium. This means that we identify po with 
p eq , which is assumed to fulfill the following equation: 



FoPeq = 0. 

The formal solution to Eq.(p 



(40) 



is 
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pi(£,t) = K [ dt'exp[T (t - t')]r lPo (Z,t') = 
Jo 

= K ( dt'exp{rot')r lPo {t, t-t') = 



K I dt' exp{Y G t')Tip eq . 
Jo 



(41) 



We assume that at equilibrium the mean value of £ van- 
ishes. As a consequence, Eq.(|4l|) yields 



where 



(t(t))=K C(t')dt', 
Jo 



C(t) = (^ea;p(rot)ri) e 



(42) 



(43) 



This concise expression means a trace over £ of the distri- 
bution obtained by applying to p eq first the perturbation 
operator T\ 1 then the unperturbed time evolution oper- 
ator exp(T Q t) and finally the variable £ itself. 

In the case of ordinary statistical mechanics one usu- 
ally makes the assumption that the variable £ has Gaus- 
sian statistics. This means a continuous variable driven 
by both the unperturbed stochastic environment and 
by the perturbation of intensity K . To define the ac- 
tion of the perturbation on the time evolution of the £- 
trajectory, let us switch off the influence of the stochas- 
tic environment. In this case the perturbation action 
would be expressed by d£/dt — K. Moving from the 
Heisenberg-like to the Schrodinger-like representation, 
we would obtain ^■ / o(£,t) = K-ggp(£,t) which implies 
that the perturbation operator is 



Ti = 



d_ 



(44) 



Note that in the Gaussian case p e g(£) is a Gaussian 
distribution proportional to exp(—^/(2(^ 2 ) eq )). Conse- 
quently C (t) of Eq. (^|) becomes 



Note that 



C(t) = 



d(x(t))/dt = (at)) 



(45) 



(46) 



and that, in the absence of perturbations, with all the 
trajectories starting at x = 0, 

<* 2 (i)} = 2(e) eq f dt' f dt"3>&"), (47) 
Jo Jo 

We adopt the subscript to indicate that the time evo- 
lution of the second moment takes place in the absence 
of perturbations. Using Eq.@, Eq.©, and Eq.@ 
supplemented by Eq.(14q), and setting equal to zero an 
arbitrary integration constant so as to ensure a response 
proportional to the perturbation strength, we arrive at 



(x(t))=K(x 2 (t)) /(2(e 



(48) 



This equation is a generalized Einstein relation |3(|] , more 
recently discussed by Barkai and Fleurov p7| . 

It is evident that the ordinary condition, (|48|), can- 
not apply in the case here under study, corresponding to 
the conditions Eqs. ( p"5| ) and (|l6|). In fact, in this latter 
case the adoption of the ordinary linear response relation 
would yield an infinite current. Earlier investigation |3§[ | 
has established that the ordinary linear response theory, 
breaks down, and a new form of linear response to per- 
turbation shows up, which, however, cannot be expressed 
in terms of unperturbed dynamics. This means that the 
response depends on the dynamical model adopted and 
must be established through numerical integration. 

In conclusion, using the strategy of Ref. p| we find 
that in the presence of feedback the free diffusion process 



dx/dt = £ 
must be replaced with 

dx/dt = + & 



(49) 



(50) 



where we have divided the current, in accordance with 
linear response theory into the sum of two terms. The 
zero-th order term is given by the fluctuation £ assumed 
to be in the same unperturbed condition as in Eq. ((f9|) . 
The first-order term is (£(t) x ). Although the variable £ 
sojourns for long times in one of the two-velocity states, 
we assume the feedback to be so weak as to be compatible 
with £ making many jumps from one velocity state to the 
other while the value of the variable x remains essentially 
unchanged. On the other hand, setting Q 



K = -A 2 x 
where A 2 is a constant, we obtain 



dx/dt = -jx(t) +£(*), 



and 
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A 2 



X, 



the susceptibility x being defined by 



X = / C(t)dt. 
Jo 



(51) 



(52) 



(53) 



(54) 



We are now in a position to address the problem under 
discussion here via three approximation steps: 
(i) We replace Eq.© with 



dx/dt — —jx(t) + r)(t), 



(55) 



where 77 is a Levy noise |39| . This means that the variable 
77 in one single time step produces jumps proportional to 
those produced by the variable £ sojourning for a time 
t = r]/W in one of the two velocity states pS|j3 0[ ]. 
(ii) We make a numerical simulation of Eq. (|52p 
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(iii) We make a completely dynamical treatment of the 
whole fluctuation-dissipation process. 

We shall refer to approximation (i) as the stochastic 
approximation. The advantage of this approximation is 
that Eq.(^) is made equivalent to an equation studied 
years ago by West and Seshadri [^.The stochastic force 
results in a phase space operator equivalent to a frac- 
tional derivative of order a, with 



P+l. 



(56) 



Thus, it is possible to show |39[ (see also Refs. [ plpO| 
for a more recent discussion of the same problemjthat 
the Fourier transform of the distribution density o~(x,t), 
a(k, t), reads: 



**(M) = -(*l*r + 7*| 



a(k,t), 



which yields the following characteristic function 



lncr(fc, t) 



(1 



b\k\ c 



which for no dissipation has the familiar form 
limlno-(M) = ~tb\k\ a . 

7— >0 



(57) 



(58) 



(59) 



However, asymptotically, in time we obtain the equilib- 
rium condition 



cr(fc, oo ) 



exp 



b\k\ a 
cry 



-b~, I k\ ' 



(60) 



The second equality defines the equilibrium parameter 
by of the Levy characteristic function. Approximation 
(ii) is referred to as the dissipative Levy walk, and has 
already been discussed in an earlier work [ p9|]3C| ], where 
it was shown that in the limiting case of 7 — > leads to 
an equilibrium distribution equivalent to that of Eq.(60). 



IV. THE NUMERICAL SIMULATION 

An interesting result of the numerical simulation of the 
dynamical model is that the equilibrium distribution is 
not uniquely determined. The condition of weak feedback 
yields a distribution which has properties different from 
those of the equilibrium distribution stemming from the 
condition of strong feedback. 

The numerical simulation of the dynamical model can 
be done both in the case of weak and strong feedback. In- 
teresting new effects are revealed by the simulation of the 
case of strong feedback. However, these are left as a sub- 
ject for future theoretical discussions. Here we illustrate 
only the simulation results concerning the case of weak 
feedback. According to the program of Ref. |l5|, we have 
developed a theory resting on the linear response theory, 
although in the non-conventional form of Section III , and 



consequently on the assumption of a very weak feedback. 
Therefore, the simulation results here discussed refer to 
a case of feedback so weak as to ensure that the require- 
ments of Section III arc fulfilled. 




-200000 -150000 -100000 -50000 



50000 100000 150000 200000 250000 



FIG. 3. The equilibrium distribution in the Levy's basin 
of attraction. The full line histograms refer to the numerical 
simulation with /3 = 0.5, 7 = 10~ 5 , W = 1, T = 50. We set 
the bin size equal to 2000. The dashed line histograms illus- 
trate the result of the numerical simulation of the dissipative 
Levy walk of Eq. (|52|. The prediction of the stochastic ap- 
proximation, or, equivalently, of the WS statistics, is denoted 
by means of the heavily dashed line. To make the figure less 
heavy we plot only the left part of the distribution predicted 
by the WS statistics. The WS equilibrium is obtained evaluat- 
ing the inverse Fourier transform of Eq. ([jo|) with i> 7 = 417771, 
a = 1.5 and 7 = 10~ 5 . This is the value of 6 7 that according 
to the theoretical prediction of Eq.(|6l|) corresponds to the pa- 
rameters of the dynamical treatment. The insert shows, for 
clarity, the enlarged portion of the figure corresponding to 
the x-axis interval [-20000,20000]. The enlargement of the 
ordinates is done after conversion to a linear scale. 

Fig. U refers to the Levy basin of attraction. We no- 
tice that the dissipative Levy walk results in pronounced 
peaks. These peaks are produced by the fact that the tra- 
jectories cannot overcross the values x = +W/X. These 
peaks signal the region within which a good agreement 
among the fully dynamic treatment, the dissipative Levy 
walk and the stochastic approximation is expected. Note 
that, as remarked in Section HI, the stochastic approxi- 
mation is equivalent to the equilibrium of the theory of 
West and Seshadri (WS) where the parameter 6 7 of 
Eq. (^) is calculated with the following formula [[3(1 : 



W((3TWY 
b-y = — — — — — sin 



709 + 1) 



(61) 



where T and (3 are defined in Eq. (|23|) and T(-) is the 
Euler gamma function. We see that in the region enclosed 
by the peaks a very good accordance among these three 
distinct approaches is found indeed. 
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FIG. 4. The equilibrium distribution in the Gaussian 
basin of attraction. The full line histograms refer to the nu- 
merical simulation with (3 = 3, 7 = 10" 4 , W = 1, f = 50. We 
set the bin size equal to 120. The dashed line histograms refer 
to the numerical solution of the Gaussian counterpart of Eq. 
([52|). The heavily dashed line is the theoretical prediction of 
the Ornstein-Uhlenbeck process with a 2 = 7.5 • 10 5 calculated 
from Eq. (^). To make the figure less heavy we plot only 
the left portion of this theoretical prediction. 

Fig. |] illustrates the situation in the Gauss basin of 
attraction. In this case the side peaks of the dissipa- 
tive walk disappear. The counterpart of the WS statis- 
tics here is the Ornstein-Uhlenbeck statistics , already 
used in an earlier publication pp[ | to predict the form of 
equilibrium distribution in the Gauss basin of attraction. 
We use the same prediction where the variance a 2 is: 



(3W 2 f 
7(/?-l) 



(62) 



and we find that the agreement with the other two treat- 
ments is again very good. 

We are convinced that the simulation results of Fig.|| 
are a reliable numerical evidence of the fact that the dy- 
namic approach to non-canonical equilibrium yields the 
analytical form proposed 18 years ago by West and Se- 
shadri |39|, here referred to as West-Seshadri (WS) statis- 
tics. However, Tsallis statistics and WS statistics have 
in common the analytical shape of the tails, which is an 
inverse power law. This might suggest that the accor- 
dance between our simulation results and Tsallis statis- 
tics is as satisfactory as, or more satisfactory than, the 
accordance with the WS statistics. We now show that 
it is not so, and that the accordance of our results with 
Tsallis statistics is much less satisfactory than with the 
WS statistics. According to the spirit of our dynamic 
approach, the comparison should be done with the Tsal- 
lis form of equilibrium corresponding to a constraint on 
the first moment of \x\. This would produce a more sig- 
nificant departure from the WS statistics, and from our 
simulation results, the WS statistics being, in fact, a Levy 
form of equilibrium statistics. The reader can easily con- 



vince him/herself about this property observing Fig [l]. 
Thus, we decided to discuss the comparison between our 
simulation results and Tsallis statistics adopting a condi- 
tion more favorable to the Tsallis statistics, namely, the 
analytical proposal of Eq. (m) . 
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FIG. 5. Tsallis and WS statistics versus the numerical re- 
sults. The WS equilibrium is denoted by the heavily dashed 
line and, as in Fig. ^, only the left part of it is illustrated. It 
corresponds to the inverse Fourier transform of the distribu- 
tion of Eq. @ with 6 7 = 417771, a = 1.5 and 7 = 10" 5 . 
The parameter 6 7 was found as in Fig. ^. The Tsallis equilib- 
rium, corresponding to the proposal of Eq. ([?]), is illustrated 
by the dotted line. The dotted line is plotted so as to coincide 
with the heavily dashed line in the region of large distances 
(6 = 6- 10~ 8 , q — 1.8). The insert shows, for clarity, the 
enlarged portion of the figure corresponding to the a;-axis in- 
terval [—20000,20000]. The enlargement of the ordinates is 
done after conversion to a linear scale. 

In Fig.|| we establish an exact accordance between 
Tsallis and WS statistics in the region of large distances. 
We see that the WS statistical prescription yields a very 
satisfactory accordance with the results of simulation also 
in the central region, while the Tsallis non-canonical equi- 
librium does not. In Fig.^ we organize the comparison 
in such a way as produce the best fitting between Tsallis 
statistics and our simulation in the central part of the 
equilibrium distribution. We see that this has the effect 
of making the discrepancy between Tsallis statistics and 
WS statistics much worse in the region of large distances. 
We see that our simulation results are much closer to the 
WS equilibrium than to the Tsallis equilibrium. We have 
to point out, furthermore, that the region of large dis- 
tances of Fig. ^, where the Tsallis form of non-canonical 
equilibrium apparently yields a better agreement with 
our numerical results, is probably a non-stationary re- 
gion, which depends on the observation time. We do not 
have yet any theory concerning this region but the result 
of the numerical observation. The numerical simulation 
shows that this region triggers non-stationary effects, due 
to the fact that the trajectories reaching large distances 
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from the central diffusion region, illustrated in Fig. 
tend to escape forever. This escape process gives rise, in 
a time scale that becomes infinitely large upon decrease 
of the feedback strength, to a diffusion process, namely 
to a form of non-equilibrium distribution. 
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FIG. 6. Tsallis and WS statistics versus the numerical re- 
sults. The WS equilibrium is denoted by the heavily dashed 
line and, as in Fig. [| only the left part of it is illustrated. It 
corresponds to the inverse Fourier transform of the distribu- 
tion of Eq. @ with 6 7 = 417771, a = 1.5 and 7 = 10" 5 . 
The Tsallis equilibrium, corresponding to the proposal of Eq. 
(0), is illustrated by the dotted line (b = 2 ■ 1CT 8 , q = 1.8). 
The dotted line is plotted so as to get the best fitting with the 
heavily dashed line in the central region. The insert shows, 
for clarity, the enlarged portion of the figure corresponding to 
the x-axis interval [-20000, 20000]. The enlargement of the 
ordinates is done after conversion to a linear scale. 

We are convinced that the WS is closer to the results 
of numerical simulation than the Tsallis non-canonical 
equilibrium. We see, however, that in the case of Fig. § 
the departure of the Tsallis prediction from the results 
of the numerical simulation is not so marked as to rule 
out this theoretical proposal. In conclusion, if we use 
only fitting arguments, we can select the constraint on 
the second moment rather than that on the first, as the 
adoption of theoretical arguments would suggest us to 
do. Eventually, the situation would appear as favorable 
to Tsallis statistics as illustrated in Fig. @. 



V. CONCLUDING REMARKS 

The present paper is not just about obtaining the best 
fit to a numerical simulation. We aim at deriving equi- 
librium from dynamics with no use of thermodynamic 
arguments. The result of our investigation is that equi- 
librium is dictated by the WS statistics rather than by 
the the Tsallis statistics. This yields a satisfactory agree- 
ment with the results of numerical simulation with no 
fitting parameters. The advocates of Tsallis statistics 



might reach a result as satisfactory as that provided by 
the dynamic theory (even if we think that actually it is 
less satisfactory) only using fitting parameters. This pa- 
per results in two important facts. The first is a numeri- 
cal simulation yielding a non-canonical equilibrium. The 
second is a theory to account for this non-canonical equi- 
librium. The theoretical arguments lead us to conclude 
that the non-canonical equilibrium must correspond to 
the WS statistics. 

If only fitting arguments were to be used, it would be 
more difficult to make a choice between the Tsallis and 
the WS statistics. However, if we invoke also theoretical 
arguments, we find that no room is left for the Tsallis 
statistics in the case of the dynamic model studied in the 
present paper. At least within the range of this dynamic 
model, we show that the non-canonical equilibrium is 
possible in nature, but it must correspond to the theo- 
retical proposal of West and Seshadri |39| rather than to 
the predictions of non-extensive thermodynamics [Q. It 
is interesting to remark that in the limiting case of very 
weak feedback our model becomes identical to the Levy 
flights subject to the Hookean force of the recent work of 
Jespersen, Metzler and Fogedby ^] , and that both mod- 
els are equivalent to that originally studied by West and 
Seshadri p9| ]. Also these authors Q found the Tsallis 
statistics to be incompatible with non-Gibbsian nature 
of the corresponding stationary solution. 

In a lucid discussion Lebowitz |^5| has recently restated 
the point of view of Boltzmann, establishing the micro- 
scopic origin of irreversible macroscopic behavior. In his 
view the adoption of the laws of large numbers is essen- 
tial, and the role of deterministic chaos becomes impor- 
tant only if it applies to a macroscopic number of non- 
interacting particles. According to Lebowitz, mixing and 
crgodicity are notions that are "unnecessary, misguided 
and misleading." In other words, this opinion reflects 
the conviction, mirrored by the handbooks of statistical 
mechanics, that the unification of mechanics and ther- 
modynamics rests on the condition N — > 00, where N 
denotes the number of degrees of freedom of the system 
under study. These statements do not, however, conflict 
with the results of Ref. |15| reported in Section [B. In 



fact, the signature of the dynamical origin of thermody- 
namics, as expressed by Eq.(|Tl|), is lost in the limiting 
case of the thermodynamic limit N — > 00, thereby mak- 
ing the controversy between the advocates of mixing and 
the advocates of N = 00 difficult, if not impossible, to 
substantiate with experimental arguments. This means, 
in other words, that the canonical equilibrium distribu- 
tion can be derived using simple arguments, based only 
on probabilistic concepts and the law of large numbers, 
or, if we wish, also on the dynamical arguments of Ref. 
p5[ . If the criterion of simplicity is adopted, one might 
be tempted to choose the former approach which leads to 
the wanted result with little or no effort, while the second 
approach yields the same conclusion after many compli- 
cated calculations based on assumed dynamical proper- 
ties, that only in a few cases have been rigorously proved. 
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The present paper shows how to extend to the case 
of boosters with no finite time scale the program of Ref. 
p5| . In this new case the dynamical approach to equi- 
librium yields a non-canonical equilibrium which is that 
advocated many years ago by West and Seshadri ppl . As 



pointed out by the theoretical discussion of Section 
it must be remarked that Eq. (|55|), yielding the form of 
equilibrium of West and Seshadri, is the consequence of 
a form of linear response, departing from the traditional 
wisdom behind the generalized Einstein relation j36|. It 
is remarkable that the only plausible form of linear re- 
sponse, resting on dynamics, yields the same prescrip- 
tion as that suggested by a phenomenological approach. 
Thus, while we do agree with Rajagopal and Abe about 
the fact that the canonical equilibrium is not the only 
acceptable form of equilibrium, we depart from them on 
the specific form that t his e quilibrium will take, since, as 
we have seen in Section [ A , the Levy statistics must not 



be confused with the generalized canonical distribution 
of Tsallis. 
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